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Abstract 

We introduce a new class of actions for staggered quarks in lattice QCD which 
significantly reduce flavour symmetry violations in the pion mass spectrum. 
An action introduced by the MILC collaboration for the same purpose is 
seen to be a special case. We discus how such actions arise from a systematic 
attempt to reduce flavour symmetry violations in the weak coupling limit. It is 
shown that for quenched lattice QCD at 6/g^ = 5.7, representative actions of 
this class give a considerable reduction in flavour symmetry violation over the 
standard staggered action, and a significant reduction over what is achieved 
by the MILC action. 
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I. INTRODUCTION 



Lattice QCD simulations have always been limited by the requirement that the lattice 
be large compared with the correlation lengths (in particular with the pion compton wave- 
length) of the theory, while the lattice spacing be small enough for physical observables to 
exhibit the scaling properties dictated by asymptotic freedom. In addition the symmetries 
of the lattice theory should approximate the continuum Lorentz and flavour symmetries. 
With the standard action this requires a lattice with a large number of sites. For this reason 
there has been considerable effort to find improved actions which obtain the desired results 
with appreciably larger lattice spacings. 

Two methods have emerged for producing such improved actions. (For a recent 
summary of such methods and a more complete set of references see In the first method, 
the action is improved in powers of the lattice spacing a, by adding higher dimensional 
operators 0. The coefficients can be calculated perturbatively or non-perturbatively 
The second method uses renormalization group analyses to determine the fixed point, 
"perfect" action [^] (This paper references to the earlier literature). For the pure gauge 
sector of lattice QCD, actions of the Liischer-Weisz form and truncated perfect actions 
1^ work very well. For the fermion part of the action, Sheikholeslami and Wohlert have 
identified the operators which contribute at low orders in a |0. For Wilson fermions both 
methods based on expansions in a and on perfect actions have had considerable success with 
lattice spacings as large as ~ 0.5 fm ||,|]. 

For staggered fermions, the improvement of the action in powers of a was discussed 
by Naik . More recently Luo has enumerated all operators which contribute at low orders 
[TT|| . Perfect action methods have also been applied to staggered fermions ||12|,|13|. While 



m a 



simple actions of both these classes show improvements over the standard action, they do 
not significantly reduce the flavour symmetry violations in the pion spectrum |T^ , p!5| . These 
flavour symmetry violations are a major part of the reason one is forced to use small lattice 
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spacings when using staggered fermions. 

The reason why these simple improvements to the staggered quark action fail to make 
significant improvements to the flavour symmetry, is that they concentrate on making the 
fermion dispersion relations near the corners of the Brillouin zone (where the components 
of the fermion momentum are close to or vr in lattice units) more continuum-like, and the 
interactions of these fermions with low momentum gauge fields also more continuum-like. 
However, the reason for the flavour symmetry violations in the staggered fermion method 
is that the gauge fields can transfer momenta large enough to take the quark from the 
neighbourhood of one corner of the Brillouin zone to the vicinity of another, which changes 
that quark's flavour. Such large momentum transfers are not suppressed by these simple 
improved actions, and so flavour symmetry violations are not significantly suppressed. 

Recently, the MILC collaboration introduced an new staggered fermion action in 
which the gauge link is replaced by the linear combination of a single gauge link and the 
sum over the six 3-link "staples" joining the same sites []TB[. Such a replacement tends to 



smooth the interaction between the gauge and fermion fields, thus reducing the coupling 
of the high-momentum components of the gauge fields to the quarks. Not surprisingly this 
action significantly decreases the flavour symmetry violation in the pion spectrum. 

We have examined how to systematically suppress such flavour symmetry violations. 
At tree level, this leads to an action which reduces the flavour symmetry breaking by an 
extra power of . We then assume that a good choice of action beyond tree level will have 
the same form, but with different coefficients. The MILC action is then seen to be a special 
case of this more general class of actions. 

We have compared the spectrum of light hadrons, with particular emphasis on the 
pions, obtained with the standard staggered quark action, the MILC action, and a subset 
of our new actions. For a preliminary search of the parameter space for these actions 
we calculated the spectrum of light hadrons on a set of quenched gauge configurations at 
/? = 6/(7^ = 5.7 on an 8^ X 32 lattice. A summary of these results was presented at Lattice'97 



. We have since calculated the spectrum of light hadrons for the standard action, a near 



optimal MILC action and a promising choice from our new actions on quenched gauge 
configurations with [3 = Q/g"^ = 5.7 on a 16^ x 32 lattice. From these calculations we 
conclude that, for an appropriate choice of parameters, our new class of actions represents 
a significant improvement over the MILC action, and confirm that both actions represent a 
considerable improvement over the standard staggered action. 

In section 2 we describe how to systematically reduce fiavour symmetry breaking at 
tree level for lattice QCD at weak coupling, and introduce our new class of improved actions 
based on this analysis. In section 3 we present our measurements of the hadron spectrum 
on quenched configurations. Section 4 gives our summary and conclusions. 



II. IMPROVING THE STAGGERED QUARK ACTION 

In the staggered fermion transcription of quarks to the lattice, the quark field on 
site n of the lattice, ip{n), is a 3 component object — a colour triplet. It lacks Dirac or 
fiavour indices. The 4 fiavours and 4 Dirac components are associated with the 16 poles 
(per colour) of the free lattice Dirac propagator. For massless quarks, these occur when 
each component of the momentum = or vr. Since interactions, in general, change the 
momenta of the quarks, they induce mixings between the degrees of freedom associated with 
different flavours and hence break flavour symmetry. 

At tree level, one can suppress flavour mixing to higher order in a, if one suppresses 
the coupling of fermions to gluons whose momentum components are all either or vr but 
not all 0. To see how this might be done, let us for the moment ignore the requirements of 
gauge invariance. Then the quark gluon coupling term in the Lagrangian could be replaced 

by 

iip\n)rif,{n)Af,{n)ilj{n + fi) - h.c. (1) 

which again gives mixing between quark flavours. We now replace in this term by 
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^ + ^1 + D-i){2 + D2 + D_2){2 + D^ + D_s){2 + D^ + D_^)A^in) (2) 

zoo 

where 

D±,A^{n) = A^{n±u). (3) 

In momentum space this is equivalent to the substitution 

AJk) ^{l + coski){l + cosk2)il + cosk3){l + cosh)AJk). (4) 
lo 

The right hand side of this equation A^{k) as A; — > and differs from A^{k) by a factor 
of only 1 + O(a^) for \k\ = 0{a). It vanishes when any component of k equals vr and is 
suppressed by a factor of O(a^) when any component of k is within 0{a) of it. Hence 
this modification of the action would suppress the tree level flavour symmetry violations by 
0(0^), which is what we want for our improved action. 

We now return to the gauge invariant theory. The quark-gluon coupling term in the 
Lagrangian is now 

il)\n)r]^{n)U ^{n)i){n + /i) - h.c. (5) 

where = expiA^. Our ansatz for a tree-level improved action is obtained by replacing 
Afj^ by in equation multiplying out the prefactors, and replacing the products of dis- 
placement operators by appropriately symmetrized covariant displacement operators. This 
leads to the replacement f/^ —>■ where 

= ^{2 + Ef^^^ + \iD,, + D-,.)] + T.i\Dup + l{D,.p + D_,,,)] 

+ Et^^^pA + ^(^^.pA + D_^,^,)]}U, (6) 

where u, p and A are summed over ±1, ±2, ±3, ±4 with |p|, |A| all different. The 

D's are the covariant displacement operators. For example D^, is defined by 

DMM = Uu{n)U^{n + p)Ul{n + p). (7) 
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The operational definition of tlie D's is as follows. The link is displaced one unit in each 
of the subscript directions. The ends of the undisplaced link are joined to the ends of 
the displaced link by products of links over the shortest paths joining the two. We then 
symmetrize over all such paths. 

We now must check that this replacement suppresses the quark-gluon interaction 
when at least one component of the gluon momentum is close to vr and the rest are near 
to 0. Since we are interested primarily in what happens to leading order in a we can write 
f/^ ~ 1 + lA^ and 1 + iA^. We then evaluate when each component of the 

momentum k of is or vr. When /ci = = ^3 = ^4 = we find A^ = A^. At the other 
corners of the Brillouin zone we find that A^ is longitudinal, and hence decouples, which is 
the desired result. 

To go beyond tree level we assume that we can use a replacement of the same form 
as the tree replacement, i.e. 

= C{xq + 2?/o + Ei^i^^ + yi{Dpu + D_^^)] + Y^^^D^p + y2{D^^p + D_^^p)] 

V vp 

+ Y^Vi^vp^ + y3(^A..pA + ^-M.pA)]}f/M (8) 
vp\ 

with C = l/(xo + 2yo + 6xi + Vlyi + 12x2 + 24?/2 + 8x3 + 167/3) Here we could modify the 
definitions of the displacement operators D which include the index ±/i to use different 
weights, depending on the positions of the links in the ±/i directions. Indeed standard 
tadpole improvement [Q] would require such changes. We have chosen not to exercise this 
option in our choices of improved actions, since even without, this class of action has 6 free 
parameters. We note that the MILC action belongs to this class, being the special case 
where Xq = 1, Xi = a;, and X2 = x^, = Vo = Vi = y2 = = 0. 
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III. THE HADRON SPECTRUM WITH IMPROVED GAUGE ACTIONS. 



One of the most visible effects of flavour symmetry violation for staggered quarks is 
seen in the pion mass spectrum. Only one of the pions is a true Goldstone boson whose mass 
vanishes as the quark mass is taken to zero. The mass differences between the non-Goldstone 
pions is typically somewhat less than that between them and the Goldstone pion. For our 
measurements, we have chosen 7^2, the other local pion, as our representative non-Goldstone 
pion. For inverse lattice spacings ~ 1 GeV, i.e. for lattice spacings ~ 0.2 fm, this symmetry 
breaking is quite large. In the chiral {rriq — ^ 0) limit m'^^/m'^ ^ 0.5 rather than 0, a 50% 
effect. Since in the real world, fn^^/m^^ ^ 0.03, this is a major impediment to working at 
such lattice spacings. 

For our measurements we have chosen to work with quenched gauge fleld conflgu- 
rations at /5 = 5.7 where the inverse lattice spacing is ~ 1 GeV. To search the parameter 
space we performed spectrum calculations on 203 independent quenched conflgurations gen- 
erated using the standard (Wilson) gauge action on an 8'^ x 32 lattice. Hadron spectra were 
calculated using a single wall source on the [odd, odd, odd) sites of the flrst time slice of 
each conflguration, gauge flxed to coulomb gauge. The propagators for local hadrons were 
measured using point sinks. 

Our measurements used the standard staggered quark action, the MILC action and 
our improved action. The vr, 1x2, p and nucleon masses are given in table |. For the MILC 
action we measured the spectrum for u = 0.5, u = 1.0 and = 00 at quark mass rriq = 0.012, 
0.02 and 0.04. From table |, we note that the flavour symmetry breaking measured as 
(m^2 ~ ^V)/^^p appears smallest for the uj = 1.0 measurements, but that the difference 
between the 3 00 values is small. From this we conclude that = 1.0 is a close to optimal 
choice for f3 = 5.7. For our improved action, even with our restricted parameterization, 
there are 6 independent parameters. To limit our choices, we chose a 1 parameter subclass 
parameterized by x, for which x„ = x"' and ?/„ = x""*"^. This choice was influenced by tadpole 
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improvement. The results we present here are for the tree level coefficients (x = 1/2) with 
quark masses rriq = 0.012 and 0.02, for x = 1 with quark masses 0.012, 0.02 and 0.04 and for 
X = oo with quark masses 0.006 and 0.012. Here symmetry breaking appears to be smallest 
for X = 1. We note also that symmetry breaking for our best improved action is appreciably 
smaller than for the best MILC action. Since the p masses are very close, this does not 
appear to be simply due to differences in the perceived lattice spacing. 

Since, even at /? = 5.7, 8^ is a rather small spatial lattice, we have confirmed and 
quantified our results using a set of 158 quenched configurations on a 16^ x 32 lattice, also 
at (3 = 5.7. The larger lattice also permitted us to go to smaller quark masses. For this 
larger lattice we have calculated the light hadron spectrum with the standard staggered 
quark action, the MILC action with uj = 1 and our improved action with a; = 1. For these 
measurements we used a single source on time-slice 1 of each configuration. This source 
was a constant for all {odd, odd, odd) sites of an 8 x 8 x 8 cube and zero elsewhere, making 
it identical to the source we used on the smaller lattice, since this seemed to produce flat 
effective mass plots. Again we worked in Coulomb gauge and used point sinks. The masses 
we obtained from fits to these propagators are presented in table ^ 

We note that the Goldstone vr masses show very little finite size effect in going from 
an 8^ X 32 to a 16^ x 32 lattice. The 7^2 masses for the MILC and improved actions also 
show relatively little finite size effect while that for the standard action shows considerably 
more. However, we note that the errors for the standard action are large, and all the errors 
in these tables are purely statistical — no estimate of the systematic error associated with 
choice or appropriateness of fits is included — so that is not clear how significant this is. 
Similar comments can be made about any apparent finite size effects in the p and nucleon 
masses. We assume that for a 16'^ x 32 lattice, where the spatial box size is > 3 fm at this (3, 
the finite size effects will be relatively small. We refer the reader to recent, more extensive 
quenched spectrum calculations at P = 5.7 for serious finite size studies and masses with 
which our standard action masses can be compared, and also for chiral extrapolation studies 
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ALilUJN 


nig 






nip 


rriN 


btandard 


0.012 


0.3156(10) 


0.659(37) 


0.891(20) 


1.217(22) 


Standard 


0.020 


0.4001(9) 


0.709(21) 


0.916(12) 


1.291(15) 


O j_ 1 1 

Standard 


0.040 


0.5495(8) 


0.858(18) 


1.006(11) 


1.548(31) 


MILC5 


0.012 


0.3227(32) 


0.462(17) 


0.793(12) 




MILC5 


0.020 


0.4083(32) 


0.533(11) 


0.809(8) 


1.181(52) 


MILC5 


0.040 


0.5550(16) 


0.674(8) 


0.893(8) 


1.378(16) 


MIL CI 


0.012 


0.3267(33) 


0.453(15) 


0.787(11) 




MIL CI 


0.020 


0.4100(29) 


0.530(11) 


0.807(8) 


1.177(48) 


MILCl 


0.040 


0.5629(16) 


0.664(6) 


0.893(8) 


1.377(15) 


MILCu 


0.012 


0.3371(32) 


0.471(13) 


0.781(10) 




MILCu 


0.020 


0.4246(28) 


0.544(10) 


0.807(7) 


1.198(43) 


MILCu 


0.040 


0.5852(15) 


0.688(5) 


0.904(7) 


1.384(20) 


naive 


0.012 


0.3562(32) 


0.453(9) 


0.788(9) 


1.170(15) 


naive 


0.020 


0.4471(30) 


0.525(6) 


0.809(7) 


1.214(11) 


imp2 


0.012 


0.3868(33) 


0.456(7) 


0.786(8) 


1.177(13) 


imp2 


020 


n 4855('22') 


n 538('5') 


n 81 5('6') 


1 264('1fi') 


imp 2 


0.040 


0.6707(18) 


0.712(3) 


0.921(5) 


1.469(40) 


impu 


0.006 


0.3660(23) 


0.448(8) 


0.770(10) 


1.163(14) 


impu 


0.012 


0.5018(22) 


0.554(4) 


0.813(6) 


1.328(23) 



TABLE I. Hadron masses for various choices of the staggered quark action at (3 = 5.7 on an 
8'^ X 32 lattice. Standard is the standard staggered action, MILC5, MILCl and MILCu are the 
MILC action with to = 0.5, 1 and oo respectively and naive, imp2 and impu are our improved 
action with x = 0.5, 1 and oo respectively. 



111,111 • 

Just comparing (m^^ — m^)/m'^p indicates that both the MILC and our improved 
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ALilUJN 


rUg 






rap 


rriN 


O j_ 1 1 

btandard 


0.006 


0.2256(5) 


0.673(23) 


0.881(10) 


1.346(14) 


O j_ 1 1 

btandard 


0.012 


0.3145(4) 


0.723(14) 


0.909(7) 


1.392(10) 


O j_ 1 1 

btandard 


0.020 


0.4000(4) 


0.813(21) 


0.941(5) 


■1 A A C\ / C\\ 

1.442(9) 


btandard 


0.040 


0.5503(4) 


0.909(11) 


1.026(5) 


1.551(7) 


MIL CI 


0.006 


0.2233(6) 


0.415(12) 


0.795(23) 


1.152(34) 


MILCl 


0.012 


0.3122(6) 


0.458(6) 


0.819(13) 


1.194(16) 


MIL CI 


0.020 


0.3993(5) 


0.521(4) 


f\ C~> f f\ / c~>\ 

0.850(8) 


1.255(12) 


MILCl 


0.040 


0.5571(5) 


0.656(3) 


0.909(5) 


-1 t~\ r\ A f c~>\ 

1.394(8) 


inip2 


0.006 


2681 ('61 


377('6'l 


763('15') 


1.111(18) 


imp2 


0.012 


0.3743(6) 


0.454(4) 


0.808(8) 


1.184(11) 


imp2 


0.020 


0.4785(6) 


0.538(3) 


0.841(6) 


1.254(12) 


imp2 


0.040 


0.6690(7) 


0.713(2) 


0.938(3) 


1.410(9) 



TABLE II. Hadron masses for various choices of the staggered quark action at /3 = 5.7 on a 
16'^ X 32 lattice. The notation is as for table |. 



actions give a considerable reduction in flavour symmetry violation over the standard action. 
In addition our improved action gives improvement over that of the MILC action. To make 
this more quantitative, we have chirally extrapolated our masses to zero quark mass. We do 
this firstly because, since the relationship between the lattice quark mass and the physical 
{MS) mass is different for each action, and there is the ambiguity as to which observable 
should be used to determine which lattice quark masses correspond to one another, niq = 
is the same for each lattice action. Secondly, the physical u and d quark masses are small 
enough that the chiral limit is a good approximation to the real world. 

Since we have 4 quark masses for each action, we can use a 3-parameter fit for our 
chiral extrapolation. For the Goldstone pion we have chosen to fit to 

ml = aniq + hm^J"^ + cm^. (9) 
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For the second pion 112 we have used a 2-parameter fit 



ml^ = a + bm„. 



(10) 



Our pion masses and these fits are plotted in figure |l|. Since our Goldstone pion masses 
have been forced to zero in the rriq = limit, m^^ is a good measure of flavour symmetry 
violation. From our fits we obtain m^^ = 0.367(44) for the standard action m^^ = 0.114(8) 
for the MILC action and m^^ = 0.073(5) for our improved action. These would be valid 
measures of improvement if we consider that the true lattice spacing is that obtained from 
some pure gluonic observable such as the string tension, or the p mass from some unknown 
"correct" fermion action. 



X 7T, 



O 7T 



0.8 - 



0.6 



C\2 



0.4 ^ 



0.2 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



standard 

MILCl 

imp2 




0.0 

0.00 0.01 0.02 0.03 0.04 0.05 

FIG. 1. Pion masses squared as functions of the quark mass on a 16^ X 32 lattice at /3 = 5.7. 
The lines are the fits described in the text. 



If, however, we follow the MILC collaboration, and consider that each fermion action 
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determines its own lattice spacing, we need to extrapolate whatever hadron mass is to be 
used to determine this spacing to rrig = for each action separately. For the p and nucleon, 
we have used a simple linear extrapolation in rriq. (Even though such fits were not great, 
we were unable to find any 3 parameter fits which did significantly better.) These fits for 
our 3 actions are shown in figures |^ and |^. This gives nip = 0.859(9) and m^v = 1.330(14) 
for the standard action, rup = 0.788(14) and rriN = 1.116(20) for the MILC action and 
TUp = 0.758(10) and ttin = 1.090(17) for our improved action in the chiral limit. Hence 
if we use the p mass to set the scale, we get m\^/m?p = 0.497(59) for the standard action 
m^^/rn^p = 0.183(12) for the MILC action and m^^/m^ = 0.127(9) for our improved action. 




1.1 
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o 



standard 



0.7 - 



X MILCl 



0.6 - 



o imp2 



0.5 



0.00 



0.01 



0.02 



0.03 



0.04 



0.05 




FIG. 2. p masses as functions of ruq for a 16'^ x 32 lattice, with linear fits. 



As discussed by the MILC collaboration, one of the effects of improving the quark 
action is to reduce the lattice spacing as determined by the p mass. (This is evident in our 
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13 = 5.7 
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o standard 



X 



o 



MILCl 
imp2 



Q g I I I I I I I I I I I I I I I I I I I I I I I I L 

0.00 0.01 0.02 0.03 0.04 0.05 

FIG. 3. Nucleon masses as functions of ruq for a 16^ x 32 lattice, with linear fits. 

results quoted above.) Before we conclude how successful our program has been, we therefore 
need to know how much flavour symmetry violation would have been improved merely by 
making such decreases in lattice spacing, without changing our quark action. To do this 
one would need to know how large the standard action flavour symmetry violations would 
be at a lattice spacings where the p mass for the standard action would have the values the 
MILC p mass and our improved p mass have respectively at /3 = 5.7. In addition, we would 
need to know the size of the MILC action flavour symmetry violations at the lattice spacing 
where the MILC p mass has the value our improved p mass has at /3 = 5.7. We do not have 
this information. However, we know that the leading flavour symmetry violations should be 
0{a^) and the ratios of relevant a's are given by the inverse ratios of the corresponding p 
masses. This would give an estimate of ml^/rrip — 0.418(50) for the standard action at a 
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lattice spacing set by the MILC p mass, and m^^/rri^p = 0.387(46) at a lattice spacing set by 
our improved p mass. Similarly we estimate the flavour symmetry violation for the MILC 
action to be Tn^^/m^p = 0.169(11) at a lattice spacing set by our improved p mass. 

We now check this C(a^) dependence against published results [jT6[ at higher /? values, 
performing the required linear extrapolations as best we can. At j3 = 5.85 the chirally 
extrapolated p mass is 0.5676(42) from which we predict m'^^/m'^ = 0.217(26) compared 
with the value calculated from extrapolated 712 and p masses namely m'^^/m'^ = 0.267(6). 
At (3 = 5.95 the chirally extrapolated p mass is 0.4629(40) from which we predict fn'^^/mj^ = 
0.144(17) compared with the direct extrapolation m^^/rn^p = 0.170(13). Thus we conclude 
that our predictions from the assumed C(a^) dependence will be if anything lower than the 
actual values. 



IV. DISCUSSION AND CONCLUSIONS 



We have introduced a new class of single-link actions for staggered fermions which 
reduce the flavour symmetry violations from O(a^) to O(a^) at tree level, where a is the 
lattice spacing, by suppressing the coupling of high momentum gluons to quarks which is 
responsible for flavour mixing. On quenched configurations a.t (3 = 5.7 where ~ 1 GeV, 
the flavour symmetry violations for local pions are reduced by ~ 65 - 75%, over those of the 
standard action and by ~ 25 - 30% over those in the MILC action. For our improved action 
this means that the flavour symmetry violations at /3 = 5.7 are approximately the size of 
those for the standard action at j3 = 6.0, i.e. at approximately half the lattice spacing. Since 
all these actions are single link, inverting the Dirac operator is no more expensive than with 
the standard action. In fact, it is considerably less expensive with the improved actions, since 
they require many less conjugate gradient iterations to reach the same level of convergence. 
For example, at rriq = 0.006, our improved action required 1000-1050 conjugate gradient 
iterations compared with 1700-1750 for the MILC action and 2600-2700 for the standard 
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action. 

Of course, although flavour symmetry violations are one of the most important barri- 
ers to using staggered quarks at lattice spacings of ^ 0.1 fm, they are not the only barriers. 
Our improvements need to be combined with improvements to the free fermion dispersion 
relations and to the gauge action. In the case of the MILC action, work of this nature has 
recently been done by Orginos and Toussaint |^ , who have also included dynamical quarks. 



The Originos- Toussaint paper does, however, point out that although their actions 
improve the vr - 7r2, mass splitting, the point-split pions do not show as much improvement. 
Although we have not measured the spectrum of these point-split pions, we expect the 
improvement to be more uniform across the pion multiplet with our action than with the 
MILC action. The reason is that at tree level, our action uniformly suppresses all flavour 
mixings. On the other hand the MILC action, at tree level, can be adjusted to maximally 
suppress some flavour mixings, but it will at best only partially suppress the others. However, 
only explicit measurement will tell if our expectations are correct. 

Another aspect of flavour symmetry violation for staggered quarks is the extent to 
which they fail to obey the Atiyah-Singer index theorem. We have shown that the MILC 



action produces only limited improvement in this area It is to be hoped that our 

improved action might fare better. 

A more serious study is needed to determine the optimal parameters in our action. 
Our action is still far from ideal for treating light u and d quarks. A tadpole improved per- 
turbative calculation of the coefficients might be helpful, although perturbative calculations 
for staggered fermions have proved disappointing in the past. One might hope that it might 
be possible to reduce the flavour symmetry violations to 0{a'^) as suggested by the tree level 



calculations, within the restrictions of a single link fermion action. The analysis of Luo |Tl 
should be helpful in reducing the operators in our action to an independent set. 

Although it is obvious that we need to reduce the coupling of high momentum gluons 
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to staggered quarks, to reduce flavour symmetry violations, it is also important to reduce 
such coupling for Wilson quarks. In the case of Wilson quarks this is to decrease chiral 
symmetry violations. This has been addressed most recently by DeGrand who introduced 
smeared link fields into Wilson fermions calculations. Such smearing could also be useful in 
reducing effects of small instantons on Wilson fermions which have the potential for creating 
problems if they are to be used as the basis for domain-wall fermions 12^ . 
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